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Abstract 

We consider a two-dimensional Fermi liquid coupled to low-energy commen- 
surate spin fluctuations. At small coupling, the hole Fermi surface is large and 
centered around (vr, tt). We show that as the coupling increases, the shape of 
the Fermi surface undergoes a substantial evolution, and at strong coupling, 
the hole Fermi surface consists of small pockets centered at (ib7r/2, ib7r/2). 
At intermediate couplings, there exists a large hole Fermi surface centered 
at (vr, vr) as well as four hole pockets, but the quasiparticle residue is small 
everywhere except for the pieces of the pockets which face the P point. The 
relevance of these results to recent photoemission experiments in YBCO and 
Bi2212 systems, and their relation to the Luttinger theorem are discussed. 
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A.C: It is a great pleasure to present a paper for a Festschrift in honor 
of David Pines. David's contributions to many fields in theoretical physics 
are well recognized throughout the world, as well as his charming personal- 
ity. David, we wish you many more successful years in Many-Body Theory. 

The problem of fermions interacting with low-energy magnetic fluctuations has attracted 
a considerable interest over the past few years particularly in connection with high-Tc super- 
conductivity (Hertz 1976, loffe and Larkin 1989, Kampf and Schrieffer 1990, loffe and Kotliar 
1990, Lee and Nagaosa 1992, Milhs 1993, Monthoux and Pines 1994, Sachdev, Chubukov 
and Sokol 1995, Altshuler, loffe and Millis 1995, Sachdev and Georges 1995, Wen and Lee 
1995). The vast majority of experimental results for cuprates indicate that the supercon- 
ducting gap is highly anisotropic and its dominant part possesses dx2_y2 symmetry. The 
process which yields this kind of pairing is the exchange of nearly antiferromagnetic spin 
fluctuations. Photoemisson experiments performed at optimal doping (Liu, Veal, Paulikas, 
Downey and Shi 1992, Liu, Veal, Paulikas, Downey, Kostic, Fleshier, Welp, Olson, Wu, Arko 
and Joyce 1992, Dessau, Shen, King, Marshall, Lombardo, Dickinson, Loeser, DiCarlo, Park, 
Kapitulnik and Spicer 1993, Ma, Quitmann, Kelley, Almeras, Berger, Margaritondo and On- 
nelion 1995) show that electrons possess a large, Luttinger-type Fermi surface whose area 
scales asl + x where x is the doping concentration. For this Fermi surface, D. Pines and Ph. 
Monthoux (Monthoux et al. 1994) have performed a number of strong couphng calculations 
of Tc within the Eliashberg formalism, and found consistent with experiments. 

There are however several issues related to the magnetic mechanism for superconduc- 
tivity, which are still subjects of debates. Naively, one would expect that the magnetic 
mechanism should be the stronger, the closer the system is to the antiferromagnetic insta- 
bility. Experimentally, however, Tc has a maximum at much larger doping concentrations 
than those where long-range magnetic order has been observed. Moreover, for a system with 
long-range order, one can study the electronic excitations and effective pairing interaction 
within the SDW technique (Schrieffer, Wen and Zhang 1989, Chubukov and Frenkel 1992). 
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At large Hubbard U, the ordered state possesses two bands of quasiparticles with Z Ri 1/2 
separated by a gap, 2A fa U. Near half-fiUing, the valence band is nearly fully occupied, 
and the conduction band is empty. Accordingly, the hole Fermi surface has a small area 
which scales as x rather than 1+x. The photoemission studies of Sr2Cu02Cl2 (Wells, Shen, 
Matsuura, King, Kastner, Greven and Birgeneau 1995) have demonstrated that the valence 
band maxima are at (7r/2,7r/2) and symmetry related points. This is consistent with the 
theoretical studies which yield hole pockets at (±7r/2, ±7r/2) near half- filling. For this kind 
of Fermi surface, there are simply no low-energy quasiparticles in the directions of momen- 
tum space, where the d^^^yi gap has the largest amplitude. Moreover, the interaction vertex 
between fermions and spin fluctuations obeys the Adler principle, i.e., it vanishes when the 
magnon momentum equals (tt, tt). As a result, the full scattering potential between fermions, 
mediated by the exchange of spin fluctuations, does not have a peak at (tt, tt) and therefore 
one cannot get d^^^y-i pairing from the magnetic mechanism (Schrieffer 1994). 

In a mean-field SDW approach, the gap, 2A, between conduction and valence bands is 
proportional to the sublattice magnetization, i.e., it vanishes when long-range order disap- 
pears. However, this is actually an artifact of the mean-field treatment. The easiest way to 
see this is to switch on a small, but finite temperature and compute lowest-order tempera- 
ture corrections to the mean-field formulae. The corrections to the sublattice magnetization 
are logarithmically divergent which is an obvious consequence of the Mermin- Wagner theo- 
rem. At the same time, the divergent terms in the corrections to A are cancelled out. As a 
result, at finite T, 2 A remains approximately equal to U, while the sublattice magnetization 
vanishes. This result is expected in some sense because the existence of two narrow bands 
with quasiparticle residue Z f« 1/2, separated by an energy gap U, is similar to what one 
obtains by just considering the atomic limit of the Hubbard model near half-filling. From 
this perspective, the role of strong antiferromagnetic fluctuations is to transform the two 
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Hubbard levels into coherent bands of quasiparticles []. 

The consideration above indicates that one can obtain an "almost" SDW description 
of the electronic states, with 2A ^ U, even without long-range order. This is consistent 
with QMC studies (Bulut, Scalapino and White 1994, Preuss, Hanke and von der Linden 
1995) which show that the density of states in a paramagnetic phase near half-filling has two 
distinct peaks roughly separated by U. However, as the system moves away from half-filling, 
the spectral weight transferes from the upper band into the lower band, and near optimal 
doping there exists just one coherent band of quasiparticles. This implies that the features of 
an SDW solution appear gradually as the system moves from optimal doping to half-filling. 

In the present paper, we analyse how the system acquires the properties of the SDW solu- 
tion. We consider a simple model in which electrons are coupled to strong antiferromagnetic 
fluctuations by 

T-Cint =9 4+<7,a'?a,/3Cfc,/34 (l) 

where g is the coupling constant, and cij are the Pauli matrices. We assume that the 
weak-coupling limit describes optimal doping and model the system's evolution towards 
half-filling by increasing the strength of the coupling constant g. The partial justification 
for this assumption follows from the diagrammatic consideration of the near atomic limit of 
the Hubbard model which will be presented elsewhere. We will also assume that the original 
SU (2) symmetry is extended to SU{N) , substitute Pauli matrices by A^^ — 1 traceless SU (N) 
generators, and perform the calculations for — oo. In this limit, the mean-field SDW 
description becomes exact. Finally, we will assume that the dynamical spin susceptibility, 
x{q, uj) is peaked at Q = (tt, tt) and its low-energy part has the form Xu^ily i^m) = A~^{6'^ + 



^Without strong magnetic fluctuations, the two bands which emerge from the Hubbard levels are 
likely to be mostly incoherent, and there also appears, at finite t, a coherent band at about U /2. 
This result was found in infinite- studies of the Hubbard model (Rozenberg, Kotliar, Kalueter, 
Thomas, Rapkine, Honig and Metcalf 1995). 
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v^{q — + + 2'y\um\)- We will consider 5 as a parameter which does not vary with 
the coupling strength. As a further simplification, we will assume that this form of the 
susceptibility holds when Vs\q — Q\ and are both smaller than some cutoff scale, C <^ W, 
where W is the fermionic bandwidth. At higher energy scales, we simply set xuiQ^^) = 0. 
This sharp cutoff however is only used here to simplify the calculations - we will see that in 
order to obtain a strong Fermi surface evolution with increasing g we only have to assume 
that the spin susceptibility gets substantially reduced at the scales which are larger than 
W. The overall factor A is determined from / dqdu xuil-,^) = S'^/{N'^ — 1). 

We now proceed to the calculations of the quasiparticle self-energy. We will first present 
the results for 7 = when a full analytical calculation is possible, and then show how the 
results are modified due to finite damping. 

The second-order self-energy contains the Green function of an intermediate state with 
momentum k + q. Since C <^ W, we can expand the fermionic energy ek+q as et+q ~ 
ek+Q + Vk+Qq, where q = q ~ Q ■ Within this approximation, we obtained the analytical 
expression for S(A;,i7) for arbitrary ratio of = \vk^Q\/vs- It turns out, however, that the 
mechanism of the Fermi surface evolution is virtually insensitive to the value and momentum 
dependence of a^. The results for arbitrary will be presented in a separate publication, 
here we just set = 1 in which case the expression for the self-energy has the simple form 

^{k,n^) = -gl ^— (2) 

for el^Q + ^1,>C^, and 



S(A;, = -gl —— ^ (3) 



^In performing the integration, we actually set Vk+Q = '>^k+Q^ integrated over ijy from —00 to +00 
and over the circle v1_^_Qq'^ + uj'^ < C^. However, we also checked that the results change very little 
if we also restrict the integration over qy to \vk+Qqy\ < C. 
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for e|_j_Q + fi^ < C^. Here we defined = — yU, and = {2/Ny/^g\S\. The location of 
the Fermi surface follows from G~^{k, 0) = — (e^ + T,{k, 0)) = 0. 

The relevant issue for our consideration is how the chemical potential fi changes with g^. 
The value of fi is obtained from the condition that 2 / dkdu G{k, u) = N/V. We solved this 
equation numerically and present the result in Fig. 1. We found that /i remains nearly the 
same as for free particles roughly up to ge = gs , when pockets are formed (see below), and 
then starts growing and approaches ge at large couplings. 

It is instructive to consider the large ge limit in more detail. As \fi\ is large in this 
limit, the condition fi^ + e^^g > is satisfied for all momenta and Matsubara frequencies 
(though not for all real frequencies!), and the self-energy is given by (@). For the full Green 
function we then immediately obtain 

G(k Q ) = ^^"^ ~ (A) 

^ ' mn.-E,{k)){rn^-E,{k)y 

where Ei^2{k) are given by 

E.,(k)-'-'^^±{d+{^i^)f. (5) 

Performing now an analytical continuation to real frequencies, we find that the full Green 
function has two poles at = £'12 (A;) ~ \fi\±ge- It is easy to check that for both solutions the 
condition for using (|^) for the self-energy (which, for real frequencies, is, {Q"^ — ef.+Q| > C^) 
is satisfied. These two solutions obey Ei^2{k) = Ei 2{k + Q) and have exactly the same form 
as the mean-field solutions of the ordered SDW state with the effective coupling g^ playing 
the role of the SDW gap, A. This exact equivalence with the SDW formulae is indeed the 
result of restricting with only the leading terms in the expansion in C/(7e - in the absence of 
a broken symmetry, k and k + Q are not equivalent points in the Brillouin zone, and in this 
sense, the quasiparticle bands at Ei^2 are only the precursors of the conduction and valence 
bands of the ordered state. Further, it follows from (^ that the quasiparticle residue for 
each of these two solutions is Z = 1/2 + 0{t/ge). Clearly then, the hole Fermi surface is 
small and encloses the area which scales as x rather than 1 + x as was the case for small ge- 



We now present our results for the Fermi surface evolution with varying (Fig. 2). 
We found that the evolution begins with the "hot spots", where = e^+Q = 0. It is clear 
from that as S(efc+Q = 0,fi = 0) = 0, the location of the "hot spots" is not shifted by 
self-energy corrections. However, at Qe = g^^^ = {25{^/C^~+~5^ — 5))^^"^ , the Fermi velocity at 
these points turns to zero. At ge > g^\ each of the "hot spots" is split into three: one is still 
at tk = 0, while the other two are located at = ±{g^ / {y/CP~T~S^ — S)){1 — (fi'^^Vfl'e)^)^^^ at 
ge < C, and at = ±ge when ge > g'f^ = C (Fig. 2b). Besides, at ge > g^e\ there appear 
(for flfe = 1) nested pieces of the Fermi surface located between points D and D' in Fig. 2c. 
For ge only slightly larger than g^\ both D and D' are located near the magnetic Brillouin 
zone boundary. As ge increases further, these points move apart and approach the Brillouin 
zone diagonal. Finally, at ge = # = C{l + 8t{\fx\-C)/C/{2t + ^Af^ - A\t'\{\fx\ - C)), the D 
and D' points from neighboring "hot spots" merge and the system undergoes a topological 
transition when the singly-connected hole Fermi surface splits into hole pockets centered at 
(±7r/2, ±7r/2) and the large hole Fermi surface (Fig. 2d). As ge increases further, the large 
Fermi surface shrinks (Fig. 2e) and eventually disappears, and at even larger ge, the Fermi 
surface consists of just four hole pockets (Fig. 2f). 

We further discuss how the quasiparticle residue along the Fermi surface changes with 
the coupling strength. At small ge, the residue is nearly A;— independent and is close to 
Z = 1. At large ge, Z along the hole pockets is again fc— independent and is approximately 
equal to 1/2. At intermediate ge, however, we found a strong momentum dependence of Z. 
In Fig. 2e we present our results for Z when large and small Fermi surfaces coexist. We 
see that the quasiparticle residue along the large Fermi surface is very small which makes 
its experimental observation problematic. The residue along the hole pocket is also very 
anisotropic and is relatively close to 1 only in a momentum region which faces the F point. 
This is consistent with the experimental observations by R. Liu et al. who reported a small 
Fermi surface in YBa2Cu30(i,s, but could only detect the Fermi surface crossing between 
the r and (7r/2,7r/2) points (Liu et al. 1992). 

We now discuss how the results above change when we include the damping of spin 
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fluctuations. We found tliat tlie general scenario of the Fermi surface evolution does not 
change, but the value of g[^^ where each "hot spots" splits into three is in fact substantially 
larger than in the absence of damping. Specifically, we found (still, assuming for simplicity 
that ttk = 1) that gP{j) = gi^\j = 0)^(7/5) where ^(x) = 1 +a;/(37r) + 0{x^) for a; < 1, 
and \Ef(x) = (16 logx/(7rx))^^/^ for x <^ 1. Near optimal doping (which corresponds to 
weak/intermediate couplings in our consideration), spin fluctuations are overdamped, i.e., 
7^5. In this case, g^^^ oc (7C/ log 7/5)^''^ which is substantially larger than g^^^ oc {5Cy^'^ 
which we obtained in the absence of damping. Though the actual numbers depend on the 
details of the band structure and specific predictions about Vs, etc, it is essential that this 
critical value does not contain 5^/^ as an overall factor. 

Now about vertex corrections. We performed the same lowest-order computations as for 
the self-energy. At small and moderate ge, the result strongly depends on the momenta of 
the external particles. In particular, it has been shown (Altshuler et al. 1995, Chubukov 
1995) that if we set the momentum in the spin propagator q = Q, and choose the incoming 
and outgoing fermions to be right at the original "hot spots" (such that both fermions are 
simultaneously at the Fermi surface), then vertex corrections in fact increase the pairing 
potential and favor (i— wave pairing. However, if we choose fermions to be at the new "hot 
spots" which appear at ge > gi^\ then the vertex corrections have opposite sign and decrease 
the pairing potential. At large ge, only a single new "hot spot" survives. In this limit, we 
found (for = 1 and to leading order in C/ge) that the full vertex function is 

nfuiifi, ^ _ ^ {inm - Ei{k)){inm - E2{k)) 

ge {K,\l,.m) — ge TTTZ - w = ^ 1,0 j 

(,^"m — ^kjKt'^l'm — ^k+Q) 

where £'1,2 are given by (^. This form of the vertex function is the same as in the mean-field 
SDW solution for the ordered state. An analytical continuation of (^) to real frequencies 
requires care as (||) is valid only when -el\> C"^ and - el_^_Q\ > C^. At Q = Ei_2(fc), 
these two conditions are satisfied, and we see that the full vertex function vanishes at 
resonance. Indeed, in the absence of the symmetry breaking, vanishing of is again 
a result of neglecting the subleading terms. Nevertheless, we see that at ge ^ C, vertex 
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corrections substantially reduce the effective pairing potential between fermions mediated 
by the exchange of spin fluctuations. 

Now consider higher-order diagrams. Here the key point is that at large ge, our results 
are the same as in the SDW description. Meanwhile, we found that the full mean-fleld 
SDW solution is given by the second-order self-energy diagram. Higher-order terms do not 
contribute because of a precise cancellation (at n — > oo) between self-energy and vertex 
corrections to the second-order self-energy term. In our case SDW solutions are reproduced 
only to lowest-order in C/ge, i.e. there is no complete cancellation. However, the higher- 
order self-energy terms are still small, at ge ^ C, compared to the second-order self-energy 
and can be neglected. At intermediate ge ~ C, the restriction with only the second-order 
diagram is, strictly speaking, unjustifled, but we expect that the scenario for the Fermi 
surface evolution will remain at least qualitatively the same as the one described above. 

Finally, some considerations on the Luttinger theorem. It is clearly violated at large ge- 
The source of this violation will be discussed in a separate publication, here we only notice 
that the key point of the Luttinger theorem is the proof, order by order in perturbation 
theory, that I = J dkdfl G{k,n) dJ:{k,n)/dn = 0. Meanwhile, for the SDW state, I is 
actually flnite, though a formal application of Luttinger's reasoning still yields / = to all 
orders in perturbation theory. The point is that the nonzero contribution to / comes solely 
from the region of (/c, Q) where the Luttinger expansion parameter A = A^|Go(^, fl)Go{k + 
Q, Q)| > 1, and perturbation series do not converge. We computed / and also the area of 
the Fermi surface with our self-energy. The results for the Fermi surface area arc presented 
in Fig. lb. Though at intermediate g, higher-order diagrams cannot be neglected and we 
cannot exactly reproduce the Luttinger theorem even when it should be valid, our numerical 
data indicate that the area stays nearly the same as for free fermions roughly up to g^'^^ when 
the hole pockets get separated from the rest of the Fermi surface. We beheve, though at the 
moment have no proof, that the Luttinger theorem ceases to work above g^\ 

It is our pleasure to thank S. Chakravarty, E. Dagotto, A. Finkelstein, D. Frenkel, L. 
loffe, R. Joynt, A. Kampf, R. Laughhn, A. Moreo, R Lee, A. Millis, M. OneUion, D. Pines, 
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FIGURES 

FIG. 1. (a) The chemical potential as a function of the coupling, ge- We used (in units of t), 
6 = 0.1, C = 0.3, t' = —0.45, X = 0.1. For free fcrmions, /x = /Lto ~ —1.16. The arrow indicates the 

(3) 

value of ge when hole pockets are formed, (b) The area of the occupied states vs ge- The dashed 
line is the area for free fermions. 

FIG. 2. (a) -(f) Fermi surface evolution with increasing ge- The parameters are the same as 
in Fig. 1. For these parameters, gi^^ = 0.21, gi^^ = 0.3, and gi^^ 0.82. The figures are for 

(2) (1) (3) (2) (3) (3) (3) 

9e = ge = 9e > ge , ge > ge > ge , ge > ge , ge > ge , and ge > ge , respectively. In (c), 
the nested pieces of the Fermi surface are shown in bold. In Fig. (e), we also presented the values 
of the quasiparticle residue along the Fermi surface. 
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